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Abstract 

A wetting transition occurs when the contact angle of a liquid drop on a surface changes from a 
nonzero value to zero. Such a transition has never been observed for water on any solid surface. 
This paper discusses the value of the temperature T w at which the transition should occur for water 
on graphite. A simple model, previously used for nonpolar fluids, predicts the value of T w as a 
function of the well-depth D of the adsorption potential. While D is not well known for the case of 
water /graphite, the model implies that T w is likely to fall in the range 350 to 500 K. Experimental 
search for this transition is warranted. Water wetting transition temperatures on other surfaces 
are also discussed. 
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I. INTRODUCTION 



Fundamental discoveries related to the thermodynamic laws describing adsorbed films on 
planar surfaces are still being made, despite the fact that these laws have been known for 
more than a century. An important example involves the contact angle 9 of a liquid drop 
on a surface and its relevance to wetting transitions . Young's equation, 



O-gs = CTls + VlgCOs(9), 



implies that 9 vanishes when a particular equality is satisfied: 



Cgs = &ls + (Jig. 



(1) 



(2) 



Here, a gs , ai s and ai g are the gas-solid, liquid-solid, and liquid-gas interfacial tensions, respec- 
tively. This description in terms of macroscopic surface tensions applies only at saturated 
vaporpressure (svp). A wetting transition occurs when 9 changes from a nonzero value to 
zero 00]. Analyses describing this transition were presented some 25 years ago in comple- 
mentary studies of Cahn and Ebner and Saam j?[ Q . Their research found that if a fluid 
does not wet a particular surface at low temperature, then a wetting transition is expected 
to occur at some higher temperature T w that is less than the bulk critical temperature T c ; T w 
thus separates nonwetting from wetting regimes. A drying transition (9 = tt), hypothesized 
in the early work, is believed now to be ruled out in the prevalent case of long range van 
der Waals forces 0, IllJ ■ 

Such wetting transitions were first observed some ten years ago in experiments invol ving 

11. 12. 



simple gases (He, Ne and H 2 ) on alkali metal surfaces 



soon afterwards, transitions were found with Hg on the surfaces of Ta and sapphi 
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ire 3, H 



The common origin of these phenomena is that the gas-surface interaction 



V(r) is only weakly attractive (implyin g nonwetting behavior at temperature T close to the 



triple temperature T 3 ) \U 



27L l28| . The dividing line between wetting and nonwetting 



behaviors at T3 corresponds approximately to the situation when the adsorption well-depth 
D = 3.5e where e is the well-depth of the adsorbate's intermolecular interaction; for a 
monatomic, classical fluid, this criterion is equivalent to the condition D/(ksT) ~ 6 at the 
triple point 30]. The physics of this phenomenon is qualitatively simple: at pressures close to 
svp the gas is on the verge of condensing to form a bulk liquid phase. If the surface provides 
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a very strongly attractive interaction, it nucleates the film, which grows continuously on 
the surface as the pressure P approaches its value Pq at svp. If, instead, the interaction is 
only weakly attractive, no such nucleation occurs below svp. In that case, a fluid droplet 
on the surface at svp will bead up instead of spreading across the surface. The difference 
between the two kinds of behavior is thus associated with the relative strength of cohesive 
and adhesive interactions. 

Interestingly, no wetting transition involving water has ever been observed, to the best 



of our knowledge [3JJ, |32| . This situation is somewhat surprising in view of the fact that 
water is a much-studied fluid that does not wet many surfaces; transitions are expected to 
occur for each of these cases. This paper is concerned primarily with the wetting behavior 
of water on graphite. There exist several reasons for this choice of substrate. One is that 
flat surfaces of graphite, with few impurities, are readily available in a number of forms. A 
second is that a large database exists for other films on this surface 0,0, 35]. In addition, 
the problem of water adsorption on graphite is related to that of water on and inside carbon 
nanotubes, a subject of considerable current interest 0000. Finally, and most 
importantly, water is known to not wet the surface of graphite at room temperature; different 
values of the contact angle have been reported, presumably due to surface imperfections 



HAH, 



43 



44J . Thus, the predicted value of T w is an open question. 



The outline of this paper is the following: The next section describes a simple model of 
the wetting transition, used previously with much success in predicting transition behavior 
of both classical and quantum fluids. Section IIHI discusses calculations of the potential and 
the resulting values of T w for water/graphite. Analogous values of T w for water on other 
surfaces are also discussed. Section HVl summarizes our results and conclusions. 



II. SIMPLE MODEL OF THE TRANSITION 



In a 1991 paper, Cheng et al. 



presented a so-called "simple model" that derives pre- 
dictions of T w from the adsorption potential and the T-dependent values of the adsorbate's 
density p and oi g . It achieves this goal by making a drastic assumption about the difference 
between two of the tensions appearing in Young's equation: 

o7s - cr gs = (Ti g + p J dzV(z) (3) 
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The logic of this relation is that ai s includes two contributions: the free energy associated 
with terminating the liquid (roughly equal to ai g ) and the liquid- surface interaction energy. 
In eqOl the domain of integration extends between the minimum in the adsorption potential 
(at z = z m i n ) and infinity. The qualitative justification for this relation includes some 
implicit assumptions: (a) that a very low density film is present on the surface just below 
the transition pressure, so that a gs is essentially the surface tension of the bare substrate, 
(b) that when an infinitely thick liquid is in contact with the substrate its excess free energy 
differs from that of the bulk liquid-vapor interface by the integrated potential energy, (c) 
that the function V(r) is assumed to depend on just the normal coordinate z and (d) that 
the integrated potential energy is adequately described by a "sharp kink" approximation, 
in which the fluid density rises from zero to its bulk value at position z = z min . If eq|3]is 
valid, then the wetting transition relation, eqEl can be manipulated to yield the following 
equation: 

' dzV(z) = I. (4) 



P 



r 



Equations El and |U provide an implicit prediction for T w , the temperature for which the left 
side equals the integral J, an explicit function of the adsorption system. The accuracy of 
eq E| can be tested by experiments, if the potential is accurately known, or by simulations, 
if it is not, by using a hypothetical potential. An extensive set of classical simulations were 
undertaken by Curtarolo et al., who tested the model by comparing its predictions for wetting 
at the triple point with simulation data for a wide variety of adsorption potentials j^j . The 
results revealed that the model works well, overall, except for situations involving the very 
least attractive interactions, in which case the error was not ne gligi ble. More recently, path 
integral Monte Carlo simulations were undertaken by Shi et al. |46j. They found the simple 
model's predictions of T w to differ from the exact results by about 10% (near 20 K) for the 
isotopes of hydrogen on several alkali metal surfaces. Comparison with experimental data 
indicates in most cases that eq El pre dicts T w relatively well when the potential is accurately 



;es in most cases that eq E Q pre dict; 
H 12, 13, 14, 15, 



known j 

One model potential that has been used extensively to treat the wetting problem is: 

This 3-9 potential is analogous to the Lennard- Jones 6-12 interatomic potential and is usu- 
ally adopted for similar reasons; it combines a rigorously correct form of the asymptotic 



attraction (V ~ C^/z 3 ) and a simple power law repulsion to yield a qualitatively plausible 
and mathematically convenient functional form. The asymptotic van der Waals coefficient 
C 3 is discussed in some detail in the next section. If one inserts eq into eq El with 

Zmin = [2 C3/ (3D)] 1 / 3 , one obtains 

11 /3\ 2 ' 3 ,„ ^V3 



When combined with eqEJ this equation yields a simple relation that we employ to determine 
T ■ 

(C 3 ^ 2 ) 1/3 = 3.33^ . (7) 

We emphasize that this relation is derived by combining the simple model of the wetting 
transition with the simple 3-9 model potential. In spite of these approximations, the result 
appears to describe most of the wetting studies carried out so far, at least semiquantitatively. 
In the next section, we evaluate the predictions of the simple model for both the 3-9 potential 
and other water/graphite potentials. 



III. THE ADSORPTION POTENTIAL 



A. Dispersion interaction 

The first aspect of the adsorption potential considered here is the van der Waals dispersion 



coefficient C3. While values of C3 exist for many gas/surface combinations 47(, no values 
have been calculated previously for the particular case of water /graphite. In the case of a 
molecule, there exist two independent contributions to this coefficient: 

C3 = Cperm + Cvdw (8) 

The first term arises from the interaction between the molecule's dipole moment p and its 
image in the dielectric substrate, while the second term is due to the cou pled dipolar charge 
fluctuations of the electrons on the molecule and those of the substrate 33]. In computing 
both of these terms, we assume a Drude model form for the surface dielectric response 
function: 

e{iE) - 1 g 



s{iE) + 1 1 + (E/E s 



\2' 



(9) 
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Here, e(iE) is the dielectric function at imaginary energy iE and E s is a characteristic elec- 
tronic energy of the solid. This form is exact, with g = 1 and E s equal to the surface plasmon 
energy, for a perfect conductor. It is also exact for a low density solid comprised of Drude 
oscillators, representing the individual molecules. For many other materials, this param- 
eterization has been shown to work well in determining the various dispersion coefficients 



49, 501 to the case of a molecule 



In evaluating C pe rm, we adapt the derivation of London 
having a permanent dipole moment p oriented at an angle relative to the surface normal, 
so that p z = p cos0, taking into account the dielectric screening at low energy (relevant to 
the slow rotation of p): 

<W0 = 9p2(1 + c ° s20) - < 10 > 

In the case of a randomly oriented, or freely rotating molecule, one obtains 
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Cperm = (Cp erm (0)) = , (11) 

since the average of cos 2 <t> for a randomly oriented dipole is 1/3. In the case of wa- 
ter/graphite, g = 0.619 |3J| and p = 1.85 Debye [5Jj, so that (in this random case) C pe rm 
=104 meV-A 3 

As to the van der Waals part of the interaction, an analogous modeling of the dynamical 
polarizability of water (with a characteristic energy E water ) yields an expression 

gaE 8 

Cvdw " 8(1 + ^-)' { 1 

v I I i^water ' 

Using data from Appendix E of Bruch et al. [33], with E WSLter =18.1 eV and a rotationally 
averaged polarizability a— 1.4 A 3 , we obtain a result C v d w = 971 meV-A 3 . From eqlHJ the 
total C3 = 1075 meV-A 3 . Note that the permanent moment's contribution to this coefficient 
is about 10%. 

Using this value of C3 in eq[7[ we can predict T w as a function of the well depth, D. In 
Figure ^ we plot T w as a function of D. The surface tension and density values for water 
used to construct Figure Q were taken from Ref. 5^ . 



B. Well-depth calculations 



There have been numerous theoretical and semiempirical studies of the water-gra 



interaction, as have recently been discussed by Werder et al. and Pertsin and Grunze 



:raphite 



The values emerging from these studies vary between 7.1 and 24.3 kJ/mole. This range 
reflects the wide variety of methods and assumptions employed in these calculations. Figure^ 
indicates the corresponding values of the wetting temperature derived from some of these 
well-depths. 



C. Results for a new potential 



Each carbon atom in g raphite has a local quadrupole moment due to the symmetry of the 



crystal (55, 



56 



57 



58 



59J. Vernov and Steele developed an angle-explicit, position-dependent 
potential that accounts for the interactions of a polar molecule with the quadrupoles on the 
carbon atoms of graphite [57] . They found that the corrugation of the water-graphite po- 
tential energy surface in the x, y plane was dramatically enhanced by the inclusion of the 
dipole-quadrupole term. Zhao and Johnson have recently developed approximate expres- 
sions for the dipole-induced dipole, dipole-quadrupole, and quadrupole-quadrupole potential 
terms for polar fluid molecules interacting with graphite 60]. The expressions are angle av- 
eraged and integrated, so that they depend only on the normal coordinate, z. For water on 
graphite only the dipole-induced dipole and dipole-quadrupole terms are important. These 
are given by 



V, 



polar 



Ap s Tcp 2 

(47T50) 2 



y 



+ 



1 



3A(z + Ay 



+ 



3k B T 



1 



+ 



1 



5A(z + A)' 



(13) 



where A is the distance between the graphene layers, p s is the density of atoms in graphite, 
Eo is the vacuum permittivity, p is the magnitude of the dipole moment of water, a c is the 
angle averaged polarizability of a carbon atom in graphite, and G c is the quadrupole moment 
on each carbon atom in graphite. The polar terms in eq^Jare added to the "Steele model" 
10-4-3 potential 0, |^ to obtain the full potential. 

The numerical value of the potential depends on the specific fluid potential used for water. 
In this work we employ the Lennard- Jones parameters from the TIP4P potential j3] for 
calculating the 10-4-3 van der Waals solid-fluid contribution. For the polar terms, we use a 
point dipole having the experimental gas phase value of 1.85 Debye. The potential depends 
on temperature, due to Boltzmann angle averaging of the dipole-quadrupole term. However, 
the temperature dependence is relatively weak. The water-graphite potential at 474 K and 
a fit of this potential to the 3-9 functional form, eqEl are plotted in Figure El The the main 
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difference between the Zhao- Johnson (ZJ) potential and the 3-9 potential developed above 
is that the ZJ potential accounts for dipole-quadrupole interactions but calculation of the C3 
term does not. The dipole-quadrupole term accounts for 12% of the energy at the potential 
minimum at 474 K, whereas the induction term is 22%. Note that the 3-9 functional form 
does not accurately describe the ZJ potential. The 3-9 potential is not as deep as the full 
potential, but it has a longer range (decays more slowly). 

With this new potential, one can evaluate the simple model to predict T w in three ways. 
One is by directly integrating V(z), including its weak T-dependence, and using eqHJ A sec- 
ond is to employ the simplified expression in eq0 using a 3-9 potential fit to the theoretical 
potential; the fit values are D = 89.8 meV and C 3 = 4751 meV-A 3 . A third way to predict 
T w is to use the theoretical well-depth (9.74 kJ/mol or 101 meV) and the theoretical value 
of C3 , i.e., use the curve in Figured] The results in these three cases are T w = 474, 416, 
and 504 K, respectively. Of these three, the value T w = 474 K, obtained from the direct 
integration, is the recommended choice since it is based on the theoretical potential, without 
recourse to a preconceived model potential. The 3-9 fit and resulting low value of T w are 
particularly suspect because the fitted value of C3 is a factor of four greater than the value 
(C3 = 1075 meV-A 3 ) derived in the previous section, which has an estimated uncertainty 
of 20 % (based on previous studies of C3 values |47(). The third value coincides with the 
one predicted from the well depth calculated by Gordillo et al. 



D. Wetting on other surfaces 

Similar calculations of C3 have been carried out for water on other surfaces, with results 
presented in Table 1. For any one of these surfaces, say x, one can compute the dependence 
D X (T W ) of the well depth on T w from eq[7| by scaling from the functional dependence D g (T w ) 
for graphite, shown in Figure Q 

D x {T w ) = (^j' 2 D g {T w ). (14) 

Here, C x /C g is the ratio of van der Waals coefficients for the surface of interest to that on 
graphite. Figure Q shows this behavior for these various surfaces. 
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IV. SUMMARY AND CONCLUSIONS 



Although water does not wet many surfaces, no wetting transition for water has been 
observed. Quite general arguments predict such a transition, unless the adsorption potential 
is extremely weak (as occurs for Ne/Cs [l(3|)- In this paper we have predicted the wetting 
temperature as a function of D, using a simple model that has been found reliable for 
spherical adsorbates j^J. One reason for concern about the reliability of the prediction is 
that the method has not yet been tested by comparing with results for water per se. Water 
is greatly affected by electrostatic forces absent from the nonpolar systems. However, the 
simple physical picture underlying the model and its previous success both suggest that it 
should be quite useful in motivating and guiding experimental searches for such a transition. 

Evidently, this subject can benefit from further research in three related directions. One 
is to carry out simulations of water /graphite based on reliable intermolecular interactions 
in order to determine the accuracy of the simple model employed here. A second is to 
further refine the water /graphite interaction potential, which is very uncertain. A third is 
an experimental search for the predicted wetting transition. 

We acknowledge very helpful discussions with Ken Jordan and Erwin Vogler. We are 
grateful for support from NSF (grants 02-08520 and 03-03916 and EEC 0085480 (JKJ)). 
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graph 


BN 


Al 


Au 


LiF 


g 


0.619 


0.38 


0.98 


0.84 


0.31 


E s (au) 


0.667 


0.71 


0.473 


0.888 


0.74 


C perm (meV-A 3 ) 


104 


64 


166 


142 


51 


C vdw (meV-A 3 ) 


971 


614 


1278 


1502 


511 


C tot (meV-A 3 ) 


1074 


678 


1444 


1644 


562 



TABLE I: van der Waals interaction coefficient and other parameters relevant to wetting of water 



on indicated surfaces. Data for g and E s are taken from Bruch et al 
are computed from equations in the text. 
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FIG. 1: Wetting temperature of water as a function of the well-depth D of the adsorption potential 
on graphite (full curve) , using eq[7|and C 3 = 1075 meV-A 3 . The square is obtained from integration 
of the ZJ potential, eq !13ja s described in the text. The open circle is based on the well-depth used 
ay Gordillo and Marti |2fj], while the filled circle uses the calculated value of Feller and Jordan 
64j |. Other curves are dependences computed for the indicated surfaces. 
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FIG. 2: Potential energy function of Zhao and Johnson 60], discussed in the text, evaluated at 
T = 474 K (full curve). The dashed curve is a 3-9 potential fit to this potential. 
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